Co-evolutionnary network approach to cultural dynamics controlled by intolerance. 



O 

> 

o 
O 



Oh' 
I 

O 

o 

O 
• I— I 

Oh. 



X 



Carlos Gracia-Lazaro,''^ Fernando Quijandrfa,^ Laura Hernandez,-^ Luis Mario Florfa,'-^ and Yamir Moreno' '*-^ 

' Institute for Biocomputation and Physics of Complex Systems (BIFI), University ofZaragoza, Zaragoza 50018, Spain 
^ Departamento de Fisica de la Materia Condensada, University of Zaragoza, Zaragoza E-50009, Spain 
^ Laboratoire de Physique Theorique et Modelisation, UMR CNRS. Universite de Cergy-Pontoise, 
2 Avenue Adolphe Chauvin, 95302, Cergy-Pontoise Cedex, France 
^ Departamento de Fisica Tedrica. University of Zaragoza, Zaragoza E-50009, Spain 
^Complex Networks and Systems Lagrange Lab, Institute for Scientific Interchange, Viale S. Severe 65, 10133 Torino, Italy 

(Dated: October 10,2011) 

Starting from Axelrod's model of cultural dissemination, we introduce a rewiring probability, enabling agents 
to cut the links with their unfriendly neighbors if their cultural similarity is below a tolerance parameter. For low 
values of tolerance, rewiring promotes the convergence to a frozen monocultural state. However, intermediate 
tolerance values prevent rewiring once the network is fragmented, resulting in a multicultural society even for 
values of initial cultural diversity in which the original Axelrod model reaches globalization. 

PACS numbers: 87.23.Ge, 89.20.-a, 89.75.Fb 



I. INTRODUCTION 

The growing interest in the interdisciplinary physics of 
complex systems, has focussed physicists' attention on agent- 
based modeling iH 01 of social dynamics, as a very attractive 
methodological framework for social sciences where concepts 
and tools from statistical physics turn out to be very appro- 
priate 111] for the analysis of the collective behaviors emerg- 
ing from the social interactions of the agents. The dynamical 
social phenomena of interest include residential segregation 
10, Et], cultural globalization opinion formation 

rumor spreading iflolfTlll and others. 

The question that motivates the formulation of Axelrod's 
model for cultural dissemination is how cultural diversity 
among groups and individuals could survive despite the ten- 
dencies to become more and more alike as a result of social in- 
teractions. The model assumes a highly non-biased scenario, 
where the culture of an agent is defined as a set of equally im- 
portant cultural features, whose particular values (traits) can 
be transmitted (by imitation) among interacting agents. It 
also assumes that the driving force of cultural dynamics is the 
"homophile satisfaction", the agents' commitment to become 
more similar to their neighbors. Moreover, the more cultural 
features an agent shares with a neighbor, the more likely the 
agent will imitate an uncommon feature's trait of the neighbor 
agent. In other words, the higher the cultural similarity, the 
higher the social influence. 

The simulations of the model dynamics show that for low 
initial cultural diversity, measured by the number q of dif- 
ferent traits for each cultural feature (see below), the system 
converges to a global cultural state, while for q above a crit- 
ical value qc the system freezes in an absorbing state where 
different cultures persist. The (non-equilibrium) phase transi- 
tion lfl2ll between globalization and multiculturalism was first 
studied for a square planar geometry f% [l3[ [Till , but soon 
other network structures of social links IUSUitII were consid- 
ered, as well as the effects of different types of noise ("cul- 
tural drift") 

external fields (modeling e.g. influential 
media, or information feedback) ll20l - l23ll . and global or local 
non-uniform couplings ll24l l25ll . 



In all those extensions of Axelrod's model mentioned in 
the above paragraph, the cultural dynamics occurs on a net- 
work of social contacts that is fixed from the outset. However, 
very often social networks are dynamical structures that con- 
tinuously reshape. A simple mechanism of network reshap- 
ing is agents' mobility, and a scenario (named the Axekod- 
Schelling model) where cultural agents placed in culturally 
dissimilar environments are allowed to move has recently 
been analyzed ll26l l27ll . In this model, new interesting fea- 
tures of cultural evolution appear depending on the values of 
a parameter, the (in-)tolerance, that controls the strength of 
agents' mobility. 

A different mechanism of network reshaping has been con- 
sidered in 1I28I I29I1 . where a cultural agent breaks its link to 
a completely dissimilar neighbor, redirecting it to a randomly 
chosen agent. At variance with the mobility scenario of the 
Axelrod-Schelling model, that limits the scope of network 
structures to clusters' configurations on the starting structure 
(square planar lattice, or others), the rewiring mechanism al- 
lows for a wider set of network structures to emerge in the 
co-evolution of culture and social ties. 

In this paper we introduce in the scenario of network 
rewiring a tolerance parameter Z controlling the likelihood 
of links rewiring, in such a way that the limit Z — 1^ re- 
covers the case analyzed in ll28[ I29I1 . where only links with 
an associated null cultural overlap are broken. Lower values 
of Z correspond to less tolerant attitudes where social links 
with progressively higher values of the cultural overlap may 
be broken with some probability that depends on these val- 
ues. The results show a counterintuitive dependence of the 
tolerance Z on the critical value q^,. On one hand, as expected 
from 1I28I I29I1 . rewiring promotes globalization for high val- 
ues of the tolerance, but on the other hand, very low values of 
Z (which enhance the rewiring probability) show the higher 
values of qc- Indeed, a non monotonous behavior is observed 
in qc{Z): Our results unambiguously show that for some in- 
termediate values of the tolerance Z, cultural globalization is 
disfavored with respect to the original Axekod's model where 
no rewiring of links is allowed. In other words, rewiring does 
not always promote globalization. On the other hand, the re- 



suiting network topology depends on q, changing from a Pois- 
son connectivity distribution P{k) to a fat tailed distribution 
for q - qc. 

II. THE MODEL 

As in Axekod's model, the culture of an agent i is a vector 
of F integer variables {(Jf{i)} (/ = 1, F), called cultural 
features, that can take on q values, af{i) = 0, 1, ...,q — 1, 
the cultural traits that the feature / can assume. The N cul- 
tural agents occupy the nodes of a network of average degree 
(k) whose links define the social contacts among them. The 
dynamics is defined, at each time step, as follows: 

• Each agent i imitates an uncommon feature's trait of a 
randomly chosen neighbor j with a probability equal to 
their cultural overlap Wij, defined as the proportion of 
common cultural features, 

where 6x,y denotes the Kronecker's delta which is 1 if 
X = y and otherwise. The whole set of N agents 
perform this step in parallel. 

• Each agent i disconnects its link with a randomly cho- 
sen neighbor agent j with probability equal to its dis- 
similarity 1 — ujij, provided the dissimilarity 1 — ixJij 
exceeds a threshold {tolerance) Z, 

l-uj^j > Z , (2) 

and rewires it randomly to other non-neighbor agent. 
The tolerance < Z < 1 is a model parameter 

First we note that the initial total number of links in the 
network is preserved in the rewiring process, so the average 
degree (fc) remains constant. However, the rewiring process 
allows for substantial modifications of the network topologi- 
cal features, e.g. connectedness, degree distribution, etc. In 
that respect, except for the limiting situation of very low ini- 
tial cultural diversity q and a very high tolerance Z (where the 
likelihood of rewiring could be very low), one should expect 
that the choices for the initial network of social ties have no 
influence in the asymptotic behavior of the dynamics. 

When the threshold tolerance Z satisfies < Z < 1, 
only those links among agents with zero cultural overlap are 
rewired, so the model becomes the one studied in lEsl I29I1 . 
On the other hand, when the tolerance takes the value Z = 
1, there is not rewiring likelihood and the original Axelrod's 
model is recovered. When Z = rewiring is always possible 
provided the cultural similarity is not complete, i.e., ujij ^ 1, 
so that it corresponds to the highest intolerance. 

The usual order parameter for Axelrod's model is 
(Smax) /N, where {S,nax) is the average (over a large number 
of different random initial conditions) of the number of agents 
sharing the most abundant (dominant) culture, and N is the 
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FIG. 1 : (color online) Order parameter as a function of the variability 
q for different values of the tolerance threshold Z. N = 900, (fc) = 
4, average over 1000 replicas. 

number of agents in the population. Large values of the or- 
der parameter characterize the globalization (cultural consen- 
sus) regime. We also compute the normalized size (Stop) /N 
of the largest network component {i.e., the largest connected 
subgraph of the network). 

III. RESULTS AND DISCUSSION 

We have studied networks of sizes N = 900, 1600; averag- 
ing over 50 - 2000 replicas. The considered cultural vectors 
have F = 10 cultural features, each one with a variability q = 
5 - 10000. We studied different values of the tolerance thresh- 
old Z G (0, 1) and different values of the average connectivity 
(fc) = 4, 10, 20, 40. Each simulation is performed for N, F, 
(fc), Z, and q fixed. For the sake of comparison with previous 
results ll28ll29tl . we will present results for (fc) — 4. 

The behavior of the order parameter for different values of 
Z is seen in Fig. [T] Like in 02811 . three different macroscopic 
phases are observed with increasing values of q, namely a 
monocultural phase, with a giant cultural cluster, a multicul- 
tural one with disconnected monocultural domains, and finally 
a multicultural phase with continuous rewiring. The nature of 
the latter phase has been successfully explained in ll28ll : At 
very large values of the initial cultural diversity q, the expected 
number of pairs of agents sharing at least one cultural trait be- 
comes smaller than the total number of links in the network, 
so that rewiring cannot stop. Here we will focus attention on 
the first two phases and the transition between them. 

In figure |2] we show the size distribution of the dominant 
culture over different realizations, measured for different val- 
ues of q, at a particular fixed value of the tolerance Z — 0.5. 
In the region of q values near the transition from globaUzation 
to multiculturalism, the distribution is double peaked, indicat- 
ing that the transition is first order, as in the original Axelrod's 
model. The transition value, qc, may be roughly estimated as 
the q value where the peaks of the size distribution are equal 
in height. The estimates of the transition points for differ- 
ent values of the tolerance Z are shown in Fig. |3] The non 
monotonous character of the graph qdZ) seen in this figure 
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FIG. 2: (color online) Histograms of Smax/N for different values of 
q, and for a fixed tolerance Z = 0.5, = 900, (k) — 4. From this 
figure one gets Qc ~ 20. 
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FIG. 3: Critical value of the diversity qc versus the tolerance thresh- 
old Z, obtained from the distribution of sizes of the dominant culture. 
= 900, (fc) = 4. See the text for further details. 



reveals a highly non trivial influence of the tolerance parame- 
ter on the co-evolution of cultural dynamics and the network 
of social ties. 

Let us first consider the (most tolerant) case Z = 0.9 that, 
except for the system size N, coincides exactly with the sit- 
uation considered in ll29ll . i.e., only links with zero cultural 
overlap are rewired. As discussed in ll29ll . for q values larger 
than the critical value for a fixed network ( qdZ = 1) ~ 60), 
rewiring allows redirecting links with zero overlap to agents 
with some common cultural trait (compatible agents), so rein- 
forcing the power of social influence to reach cultural global- 
ization. Once all links connect compatible agents, rewiring 
stops |l30l ]. From there on, the network structure will re- 
main fixed, and globalization will be reached with the pro- 
viso that the network has so far remained connected. This is 
the case for most realizations (for N = 900) up to values of 
q ^ 240. Increasing further the cultural diversity q, increases 
the frequency of rewiring events and slows down the finding 
of compatible agents, favoring the topological fragmentation 
into network components before rewiring stops. Under these 
conditions, the asymptotic state will consist of disconnected 
monocultural components. 

On one hand, network plasticity allows to connect com- 
patible agents, so promoting globalization; but on the other 
hand it may produce network fragmentation, so favoring mul- 
ticulturalism. What we have seen in the previous paragraph 
is that for Z = 0.9 the first effect prevails over the second 
one up to qc{Z — 0.9) ~ 240. Going from there to less 
tolerant situations (decreasing Z), increases the likelihood of 
rewiring, making easier that network fragmentation occurs be- 
fore rewiring stops. This has the effect of decreasing the criti- 
cal value qc- In fact, from Fig.[3]we see that for Z = 0.7, 0.6, 
and 0.5 multiculturalism prevails for cultural diversities where 
the original Axekod's model shows cultural globalization. In 
these cases network plasticity promotes multiculturalism in a 
very efficient way: Agents segregate from neighbors with low 
cultural similarity and form disconnected social groups where 
full local cultural consensus is easily achieved, for q values 



low enough to allow a global culture in fixed connected net- 
works. 

For very low values of the tolerance parameter, though net- 
work fragmentation occurs easily during the evolution. Fig. [3] 
shows that globalization persists up to very high values of the 
initial cultural diversity q. To explain this seemingly para- 
doxical observation, one must realize that network fragmen- 
tation is not an irreversible process, provided links connect- 
ing agents with high cultural overlap have a positive rewiring 
probability. Under these circumstances, transient connections 
among different components occur so frequently so as to make 
it possible a progressive cultural homogenization between 
components that otherwise would have separately reached dif- 
ferent local consensuses. Fig. |4] illustrates the time evolution 
for q = 100 and different values of Z. Panel (a) shows an ex- 
ample of cultural evolution where network fragmentation re- 
verts to a connected monocultural network for Z ~ 0.2. Panel 
(b), that corresponds to Z = 0.6, shows that social fragmen- 
tation persists during the whole evolution, while in panel (c), 
which corresponds to the most tolerant situation { Z ^ 0.9), 
the network remains connected all the time. 

The degree distribution of the network is Poissonian cen- 
tered about (k) for all q values, except for q ^ qc where it 
becomes fat tailed, with several lowly connected (and discon- 
nected) sites. For very high q values, in the dynamical phase, 
the network rewiring is esentially random, so Pq{k) is again 
Poisson like, centered around (fc). 

IV. SUMMARY 

In this paper we have generalized the scenario for co- 
evolution of Axelrod's cultural dynamics and network of so- 
cial ties that was considered in [28, 29], by introducing a tol- 
erance parameter Z that controls the strength of network plas- 
ticity. Specifically, Z fixes the fraction of uncommon cultural 
features above which an agent breaks its tie with a neighbor 
(with probability equal to the cultural dissimilarity), so that, 
the lower the Z value, the higher the social network plasticity. 
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of the initial cultural diversity where multiculturalism would 
have been the outcome for fixed networks. On the contrary, 
for intermediate values of the tolerance, the network plasticity 
produces the fragmentation of the (artificial) society into dis- 
connected cultural groups for values of the initial cultural di- 
versity where global cultural consensus would have occurred 
in fixed networks. For very low values of the tolerance, so- 
cial fragmentation occurs during the system evolution, but the 
network plasticity is so high that it allows the final cultural ho- 
mogenization of the transient groups for very high values of 
the cultural diversity. Intermediate tolerances promote multi- 
culturalism, while both extreme intolerance and extreme toler- 
ance favor the formation of a global culture, being the former 
more efficient than the latter. 
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FIG. 4: (color online) Time evolution of mean overlap and number 
of topological clusters for diferent values of tolerance Z = 0.2(a), 
Z = 0.6(b), Z = 0.9(c). N = 900, q = 100. See the text for 
further details. 



Our results show that the network plasticity, when con- 
trolled by the tolerance parameter, has competing effects on 
the formation of a global culture. When tolerance is highest, 
network plasticity promotes cultural globalization for values 
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